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Using 214 fb^^ of data recorded by the BABAR detector at the PEPII electron-positron collider, we 
study the decay Df -^ K^K^e'^Ve. Except for a small S-wave contribution, the events with K^K~ 
masses in the range 1.01-1.03 GeV/c^ correspond to (j) mesons. For Of — > 4>e'^Ve decays, we measure 



the relative normalization of the Lorentz invariant form factors at q^ = Q,rv ~ y(0)/y4i(0) = 1.849± 
0.060 ± 0.095, r2 = A2(0)/yli (0) = 0.763 ± 0.071 ± 0.065 and the pole mass of the axial-vector form 
factors niA ~ (2.281q'^| ± 0.18) GeV/c^. Within the same K^ K~ mass range, we also measure the 
relative branching fraction B(_D+ -^ K+K'e+Ue)/B{D+ -^ K+K-tv+) = 0.558±0.007±0.016, from 
which we obtain the total branching fraction B{D+ -^ (t)e+v^) = (2.61±0.03±0.08±0.15) x 10^^ . By 
comparing this value with the predicted decay rate, we extract j4i(0) = 0.607±0.011±0.019±0.018. 
The stated uncertainties are statistical, systematic, and from external inputs. 

PACS numbers: 12.15.Hh, 12.38.Gc, 13.20.Fc, 14.40.Lb 



Charm semileptonic decays can help to validate predic- 
tions from lattice QCD through precise measurements of 
hadronic form factors. Such measurements have been 
performed by BABAR for the D^ -^ K~e^Ve decays [Ij. 
The Z)+ -^ 0e+i^e jj] channel is well suited to study form 
factors in semileptonic decays of charm mesons to a vec- 
tor particle because the </> meson is a narrow resonance 
which can be well isolated experimentally. Because of the 
higher mass of the spectator s-quark, form factor deter- 
minations for this process by lattice QCD are expected to 
be more accurate than for non-strange D mesons. How- 
ever, measurements of this decay mode are impacted by 
the lower production rate for D+ mesons and higher 
backgrounds. Form factors in Z?+ — > (pe^Ve have been 
previously studied by photoproduction experiments, at 
Fermilab 0,1,11,0, and by CLEOII at the CESR e+e" 
coUider also operating at the T{AS) [7|. In charm me- 
son semileptonic decays, a (j) meson is expected to origi- 
nate only from the D'^ . A possible contribution from the 
Cabibbo suppressed D+ -^ 0e+z^e decay, through the dd 
component of the meson [8| is neglected [91 . 

Using 214 fb^^ of data collected at the T(4S') reso- 
nance by the BaBAR detector, we measure the _D+ 



to S- and P-wave contributions, these form factors can 
be written as: 

^1 = •?^io+-^ii cosOk, Ti = —=Ti\, Tj, = —=Tz\- (1) 

The form factors Tij depend only on rr?xK ^^^ 9^j -^lo 
characterizes the S-wave contribution, whereas the Ti\ 
correspond to the (^ meson: 



(2) 



Tix = V3T:qH,{q ,m)A^{m), 



where the (j) meson decay amplitude Aiji{m) is taken to 
be a relativistic Brcit-Wigncr distribution with a mass- 
dependent width including a Blatt-Weisskopf damping 
factor [I4I. The form factors Hi 2.3 can be expressed in 
terms of the Lorentz invariant form factors V and A1.2 
[l3| , for which we assume a q^ dependence dominated by 
a single pole: 



Viq') = 



y(o) 



1 - q^/m{ 



AiMq') 



^1,2(0) 
1 — q^ l"ir?^ 



(3) 



mA and my are the pole masses, usually fixed to 
the values of corresponding resonance masses: rtiA — 



q" = 0, the ratios of the form factors V and A2 relative to 
A\ are denoted by ry and ti-, respectively. The S-wave 
contribution is parameterized assuming /q production: 



if+if-e+z^e channel decay characteristics, for events pro- 2.5 GeV/c2(~ to^^J and my = 2.1 GeV/c2(~ m^*). At 

duced in the continuum e^ e~ -^ cc. The analysis focuses 

on the (pe^Vf^ final state in the K^K^ invariant mass 

range between 1.01 and 1.03 GeV/c^. The ip resonance is 

dominant in this K^K~ invariant mass region although 

a small S-wave component is observed, for the first time, 

through its interference with the ip. 

The differential decay rate for Df -^ K^K^e^Vf. de- 
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Pkk mo. 
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(4) 



pends on five variables |1C 



''KK' 



the mass squared of 



the K^K^ system; q^, the mass squared of the e+^^e sys- 
tem; cos^e {cos Ok), where 6^ (Ok) is the angle between 
the momentum of the e+ {K~^) in the e+z/g (K^K^) rest 
frame and the momentum of the e'^Ve (K^K^) system 
in the D'^ rest frame; and x, the angle between the nor- 
mals to the planes defined in the Df rest frame by the 
K^K^ pair and the e^v^ pair. When analyzing a Dj 
candidate, the direction of the K~ is used in place of 
the K^ and x is changed to — %. The expression for the 
differential decay rate as a function of these variables is 
given in ref. llj. Neglecting contributions proportional 



where r^ is a normalization factor and pkk is the mag- 
nitude of the three-momentum of the K~^K~ system in 
the D^ rest frame. The values of the /o parameters 
(m/o , gTT , r^s ) are taken from Ref. [14 1 . 



A detailed description of the detector and the algo- 
rithms used for charged and neutral particle reconstruc- 
tion and identification is provided elsewhere 15|. Monte 



to the square of the electron mass, it depends on three 
hadronic form factors which are related to the three pos- 
sible helicity values of the hadronic current. Restricting 



Carlo (MC) samples of T(4S') decays, charm and other 
light quarks pairs from continuum events are generated 
using a GEANT4 [l6|. Quark fragmentation, in con- 
tinuum events, is described using the JETSET package 
[17|. Signal MC events are generated with seven times 
the equivalent statistics of the data, using a simple pole 
model for the form factors with niA = 2.5 GeV/c^ and 
mv — 2.1 GeV/c^. The simulation of the characteristics 
of Z?+ production is corrected to account for measured 



differences compared to data, 
simulated with PHOTOS [11 1. 

We reconstruct Df -^ K^K 
produced in e^e'^ — > cc events. 



Radiative processes are 



'e+i^e decays, for _D+ 
The hadronization of 
the cc system leads to the formation of two jets, emitted 
back-to-back in the center-of-mass (cm.) frame. The 
analysis method is similar to the one used for the decay 
D" -^ K~e^i'e P|. The only differences are that the cas- 
cade from a D* is not used to evaluate the signal, and the 
detector performance for the Dj" reconstruction is mea- 
sured using _D+ -^ (pTT^ decays rather than the cascade 
decay D*+ -^ D^tt+, D° -^ K-ti+. 

The event thrust axis is determined from all charged 
and neutral particles in the cm. system and its direc- 
tion is required to be in the range | cos(0thrust)| < 0.6 
to minimize the loss of particles in regions close to the 
beam axis. A plane perpendicular to the thrust axis is 
used to define two hemispheres, equivalent to the two jets 
produced by quark fragmentation. In each hemisphere, 
we search for the decay products of the Df , namely a 
positron, of momentum greater than 0.5 GeV/c, and two 
oppositely charged kaons. Since the v^ momentum is un- 
measured, a kinematic fit is performed, constraining the 
invariant mass of the candidate K^K'e^Ve system to 
the Dj" mass. In this fit, the D+ direction and the neu- 
trino energy are estimated from the other, charged and 
neutral, particles measured in the event. The D+ di- 
rection is taken as the direction opposite to the sum of 
the momenta of all reconstructed particles, except for the 
kaons and the positron associated with the signal candi- 
date. The neutrino energy is estimated as the difference 
between the total energy of the jet containing the can- 
didate and the sum of the energies of all reconstructed 
particles in that hemisphere. The 13+ candidate is re- 
tained if the x^ probability of the kinematic fit exceeds 
10-2. 

Sizable backgrounds arise from T(4S') -^ BB decays 
and two-jet events from e+e" -^ qq, q = u, d, s, c. 
Backgrounds are predominantly rejected by using two 
Fisher discriminant variables that exploit differences in 
the production characteristics of hadrons in signal and 
background. The first variable is used to separate signal 
in jet-like cc events from BB with a more spherical topol- 
ogy. The chosen cut retains 71% of the signal and rejects 
86% of the BB background. The second Fisher discrim- 
inant uses variables related to the different production 
characteristics of particles from Dg decays and c-quark 
fragmentation. The selected cut retains 71% of the signal 
decays, and rejects 72% of the background. The overall 
signal efficiency is approximately 4.5%. Figure[T^) shows 
the K^K^ invariant mass distribution for the selected 
decays compared to the simulation. There are 31,839 
events in the signal region, with an estimated background 
of 20.3%. About 70% of the total background is peaking, 
corresponding to a decay combined with an electron 
from another source. The interference between S- and P- 
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FIG. 1: a) K'^ K~ invariant mass distribution from data 
and simulated events. MC events have been normalized to 
the data luminosity according to the different cross sections. 
The arrows indicate the selected K'^ K~ mass interval. In b), 
each event is weighted by the measured value of cos Ok ■ Neg- 
ative entries are produced by the cc background asymmetry 
in cos 9 k ■ 



waves generates an asymmetry in the cos Ok distribution 
which is revealed in Fig. [Ud), where events have been 
weighted by cos 9k ■ 

To extract Ns (the number of reconstructed signal 
events), rv^ ?'2, friA and rg, we perform a binned maxi- 
mum likelihood fit to the four-dimensional decay distri- 
bution in the variables q^ , cosOe, cos 9k and x- The 
sensitivity to mv is weak and we fix this parameter to 
2.1 GeV/c^. The data are divided into 625 bins, with five 
equal-sized bins per variable, and 

625 

/: = -^ln7'(nf'-|nfC). (5) 

For each bin i, V {n'^^^'^\n^'^) is the Poisson probability 
to observe nf^'^'^ events when n^*-" events are expected. 



nf^iX) = Ns 



WtotW 



,BMC 



(6) 



Here rq ^ is the number of signal MC events with re- 




TABLE I: Measured values of the parameters Ns, ry, r^, 
niA (GeV/c^), ro (GeV~^), their statistical and systematic 
uncertainties. 



FIG. 2: Projected distributions of the reconstructed values of 
the four kinematic variables. The data (points with statistical 
uncertainties) are compared to the sum of four histograms, 
the fitted signal (light hatched) and the estimated background 
contributions (dark shaded). 



constructed values of the four variables corresponding to 
bin i and nf^'^ is the number of estimated background 
events. They are obtained from MC simulation, corrected 
for measured differences between data and simulation. 
Weights, Wj, are evaluated for each event, using the gen- 
erated values of the kinematic variables, thus accounting 

for resolution effects. Wtot (A) = ^j^i Wj (A) is the sum 
of the weights for all simulated signal events (TV^^^^) and 
A corresponds to the parameters to be fitted. The data 
and results of the fit are shown in Fig. [5] and listed in 
Table H] From the fit we extract a contribution due to 
S-P wave interference. The value obtained for ro corre- 
sponds to a S-wave fraction of (0.22^Q;Jg)% of the decay 
rate. 

In the fitting procedure, two sources of statistical fluc- 
tuations are not included. They originate from the finite 
sample of simulated signal events and the estimate of the 
average number of background events in each bin. These 
effects are evaluated with parameterized simulations and 
included in the systematic uncertainties. Other system- 
atic effects have been assessed to account for the un- 
certainties in the c-quark hadronization, the background 
contributions, and the remaining uncertainties in the 
simulation of the detector response. They are summa- 
rized in Table m 

Corrections to the simulation of the c-quark fragmenta- 





Ns rv r2 ruA 


To 


Fitted value 


25341 1.849 0.763 2.28 


15.1 


Statistical uncertainty 


178 0.060 0.071 +l'it 


2.6 



MC statistics 81 

cc fragmentation 5 

Background subtraction 480 

Detector effects 12 



0.029 0.034 0.09 0.0 

0.031 0.014 0.07 0.4 

0.081 0.049 0.14 0.6 

0.021 0.021 0.07 0.6 



Total syst. uncertainties 488 0.095 0.065 0.18 1.0 



tion were performed iteratively, comparing variables used 
in the event selection for samples of D'^ -^ (f)Tr^ decays 
and applying a weight which depends on the values of 
these variables. We adopt the observed changes in the fit 
parameters for the last step in this iterative process as an 
estimate of the systematic uncertainty. Furthermore we 
assume a 30% uncertainty in the simulation of radiative 
effects. 

The peaking and combinatorial background compo- 
nents from e^e~ — > cc events have been studied sepa- 
rately. The peaking background contributions are stud- 
ied by measuring inclusive (j) production in events with 
a fully reconstructed D*'^ or Df decay. The combinato- 
rial background consists mainly of events with a charged 
lepton, one kaon from a D decay and a second kaon 
from fragmentation. We have measured the rate, mo- 
mentum and angular distributions of K^ accompanying 
a D^, D*~^ or £)+ meson in data and corrected the cor- 
responding simulation. 

After these corrections, the K^K^ distribution for se- 
lected signal events in MC and data agree to within 
10% above 1.03 GeV/c^, and this remaining difference 
is adopted as the uncertainty in the normalization of the 
combinatorial background. The BB background is ob- 
tained from the difference of the data recorded at the 
T{AS) resonance and the data recorded 40 MeV/c be- 
low. The related systematic uncertainties are obtained 
from the statistical accuracy of these measurements and 
from the uncertainty (0.25%) between the relative nor- 
malization of the two data samples. Systematic uncer- 
tainties also originate from the simulation of the detector 
response. There are small differences in the efficiencies 
for charged particle reconstruction and electron and kaon 
identification. They lead to data-MC differences in the 
reconstruction of the D'^ direction and the neutrino en- 
ergy. They are estimated using D'^ — > ^tt"*" decays. 

We measure the D^ -^ K^K^e^Vf. branching frac- 
tion relative to the decay, D~^ -^ K'^K~tt^ for which we 
adopt the K+K' mass interval, 1.0095-1.0295 GeV/c^, 
to match the range used by CLEO-c for the D"^ -^ 
K^K^TT'^ branching fraction measurement [l9|. Specifi- 
cally, we compare the ratio of rates for the two channels 



in data and simulated events so that most systematic un- 
certainties cancel [1] . In the considered mass intervals, we 

obtain Rds = ^if^°V^^"^fJ"''r:^ = 0.558 ±0.007 ±0.016. 
Systematic uncertainties are summarized in Table [III 
They originate mainly from selection criteria that are 
not common for the two channels. Differences in the im- 
pact of the two Fisher discriminants have been estimated 
by varying the selection cuts and differences in particle 
identification for electrons and pions are accounted for. 
The uncertainty on Ns is taken from the previous fit; it 
is dominated by uncertainties in the background evalua- 
tion. We translate the ratio R^s to a branching fraction, 



TABLE II: Summary of the relative systematic uncertainties 
on Rds- 



Source 




Relative variation 


Fisher variable against 


cc events 


±1.77% 


Fisher variable against 


66 events 


±0.58% 


Fitted signal (Ns) 




±1.92% 


FID corrections 




±0.74% 


Df production 




±0.20% 


Mass constrained fit 




±0.61% 


Total systematic uncertainty 


±2.85% 



using B{D+ -^ K+K-TT+) = (1.99 ± 0.10 ± 0.05)% [19|, 
correcting for the finite mass range used to select signal 
events (86.37 ± 1.22)%, subtracting the S-wave contribu- 
tion, and taking B{(t> -^ K+ R-) = (49.2 ± 0.6)% 8]. 
We find: 



B{Di 



he+ve) = (2.61 ± 0.03 ± 0.08 ± 0.15) x 10" 



where the last quoted uncertainty corresponds to external 
inputs. 

In conclusion, we have studied the decay Df -^ 



K+K- 



e^Vt 



with a sample of approximately 25,000 sig- 



nal events, which greatly exceeds any previous measure- 
ment. This decay is dominated by the 4> vector meson; we 
measure a small S-wave contribution, possibly associated 
with /o -^ K+K-, corresponding to (0.22to:o8 ± 0.03)% 
of the K^K^e^Vf. decay rate. We have extracted form 
factor parameters from a fit to the four-dimensional de- 
cay distribution, assuming single pole dominance and ob- 
tain: rv = y(0)/Ai(0) = 1.849 ±0.060 ±0.095, r^ = 
^2(0)Mi(0) = 0.763 ± 0.071 ± 0.065 and the pole mass 
of the axial-vector form factors niA = (2.28lQ'^g ± 
0.18) GeV/c^. For comparison with previous measure- 
ments we also perform the fit to the data with fixed pole 
masses rn^ = 2.5 GeV/c^ and my = 2.1 GeV/c^, ignor- 
ing also the small S-wave contribution. The results on ri 
and Ty represent a large improvement in statistical and 
systematicprecision compared to earlier measurements 
BSSSIS (see Table [nH. 
We also measure the relative branching fraction 



TABLE III: Results from previous experiments and present 
measurements. 



Experiment 


rv 


r2 


E653 f3 


2.3ti-^±0.4 


2.11-^ ±0.2 


E687 [4 


1.8 ±0.9 ±0.2 


1.1 ±0.8 ±0.1 


E791 f5 


2.27 ±0.35 ±0.22 


1.57 ±0.25 ±0.19 


FOCUS [6] 


1.549 ±0.250 ±0.145 


0.713 ±0.202 ±0.266 


CLEOII [7] 


0.9 ±0.6 ±0.3 


1.4 ±0.5 ±0.3 


BABAR 


1.807 ± 0.046 ± 0.065 


0.816 ± 0.036 ± 0.030 



0.007 ± 0.016, from which we obtain the total branch- 
ing fraction B{D+ -^ cj)e+Ve) = (2.61 ± 0.03 ± 0.08 ± 
0.15) X lO^'^. By comparing this quantity with the pre- 
dicted decay rate, using the fitted parameters for the 
form factor pole ansatz we extract ^i(O) = 0.607 ± 
0.011 ± 0.019 ± 0.018. Here the third uncertainty refers 
to the combined value from external inputs, namely the 
branching fractions of the Z5+ into K^K^tt^ and of the 
(j) into K+K~, the D+ lifetime [(500 ± 7) x lO^i^s] and 
Vcs = 0.9729 ± 0.0003. Predictions for this decay chan- 
nel of lattice QCD calculations, in the quenched approx- 
imation d^, give: ry = 1.35+J]:S]^, ra = 0.98 ± 0.09, 
niA = 2.42l°;22 GeV/c^ and yli(O) = 0.63 ± 0.02. They 
agree with our determination of ^i(O), r2 and uia, but 
are lower than the measured value of ry . The measured 
form factor's ratio r2 is in agreement with the value ob- 
tained for the same parameterization for the vector de- 
cay D —^ K*e^Ve, whereas ry is two standard deviations 
higher [8| . The branching fraction presented here agrees 
well with the value (2.68 ± 0.13) %, consistent with the 
assumption of equal semileptonic decay widths for the 
different charm mesons. 
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K+K-e+i'e)/B{Df -^ K+K- 
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